American Journal of Biological Anthropology

Maya peoples have undergone several largescale demographic transformations over the past 3000 years, including the formation and depopulation of large urban centers, climatic events, and multiple colonial regimes. Archaeological and stable isotope evidence from across the Maya lowlands throughout these periods paints a picture of a demographically mobile population with complex and dynamic relationships among some of the largest urban polities in antiquity. The nature of the genomic relationship between ancient Maya urban centers, as well as their relationship with present-day Maya communities, however, remains under-characterized. Here, we present mitogenomic data of 16 ancient Maya Ancestors recovered from interments in the ancient Maya cities of Caracol and Santa Rita Corozal located in Western and Northern Belize. These ancient data were combined with 110 newly collected contemporary individuals from Northern Belize, along with nearly 600 previously published ancient and contemporary comparative mitogenomes recovered from sites spanning Mesoamerica and the Caribbean to reconstruct paleodemographic and mobility patterns of maternal lineages through time. Our results are consistent with broadly interconnected Maya populations that have maintained maternal genetic continuity across much of the Holocene, and which are not structured by linguistic family.

The relative roles of genetic drift, inbreeding, and natural selection in island evolution have been fundamental themes in evolutionary biology for nearly a century. However, empirical studies have yet to reach a consensus, and research integrating both genotypes and phenotypes in primates remains limited. The present study aims to estimate the demographic history of an island primate population, the Yakushima macaque (Macaca fuscata yakui), and evaluate the genomic and morphological consequences of its isolation. To this end, we conducted population genetic analyses based on whole-genome sequencing and geometric morphometric analyses of cranial shape. Our results indicate that M. f. yakui has been isolated with an effective population size of several thousand individuals for over 100,000 years. Although we detected signatures of purifying selection, deleterious alleles have not been sufficiently purged. The degree of fluctuating asymmetry in males was slightly, yet significantly, higher in M. f. yakui than in mainland populations, whereas this trend was not observed in females. Marginal evidence of directional selection was found for cranial size, but not for shape. Overall, our findings suggest that due to long-term isolation, genetic drift and inbreeding have exerted a substantial influence on shaping genotypic and phenotypic variation in this island population.

Introduction: Killer-cell immunoglobulin-like receptor (KIR) gene content shapes natural killer-cell responses and varies substantially among populations. We examined whether worldwide KIR diversity reflects regional structure, isolation by distance, serial founder effects, or latitudinal gradients. Data and methods: Carrier frequencies for 11 variable KIR genes were obtained from the Allele Frequency Net Database. The primary analysis included 189 locally sampled, non-diaspora, non-mixed populations from 70 countries, comprising 25,469 individuals; 161 populations had complete 11-gene profiles. Diversity indices, ordination, clustering, permutation-based variance decomposition, Mantel tests, and regression models were used to characterize geographic structure. Serial founder effects were evaluated using four distance measures from East Africa, both globally and within the Old World. Models of latitude incorporated spatial autocorrelation. Prespecified sensitivity analyses evaluated alternative population sets, gene panels, geographic classifications, weighting schemes, and multivariate outlier exclusion. Results: Geographic region explained 53.5% of variation in KIR profiles, substantially exceeding the permutation-null expectation of 6.9%, although multivariate dispersion differed among regions. Twenty-five of 28 pairwise regional comparisons remained significant after false-discovery-rate correction. Isolation by distance was modest (Mantel r=0.3003) and weakened after controlling for region (r=0.1557). KIR diversity did not decline with distance from East Africa globally but decreased strongly within the Old World (expected heterozygosity {beta}=-0.0066 and standardized gene content {beta}=-0.0789). B-content decreased nominally with absolute latitude, but latitude was collinear with region and the association survived neither spatial modelling nor sensitivity analyses. Regional differentiation remained robust across analytical specifications, whereas the Old World decline was carried by gene content once standard errors were clustered by country. Conclusion: Worldwide KIR gene-content variation shows pronounced regional differentiation, modest isolation by distance, and no independent latitudinal gradient. Diversity declines with distance from East Africa across the Old World, supporting serial founder effects, but not once the Americas and Oceania are included.

Fossil central nervous systems (CNS) have illuminated several facets of total-group euarthropod evolution, including their phylogenetic affinities, the origin of the head, and their ancestral neuroanatomical organization. However, the CNS remains unknown in trilobites, the most diverse and earliest group of Palaeozoic euarthropods, but can potentially clarify fundamental questions that remain disputed despite thousands of described species and countless studies spanning over two centuries, such as their affinities relative to modern representatives. Here, we describe the CNS of the cheirurid Anacheirurus adserai from the Lower Ordovician Fezouata Shale. Anacheirurus features a brain with well-developed proto, deuto- and tritocerebral neuromeres, an anteriomedian protocerebral protrusion, a central complex, an esophageal foramen, and a ganglionated ventral nerve cord. The CNS is replicated in framboidal iron that differs in texture and composition from the exoskeleton and rock matrix. The overall organization of Anacheirurus's CNS most closely resembles that of Cambrian stem-group mandibulates, and critically shares the presence of a central complex with a protocerebral bridge and lateral accessory lobes with extant pancrustaceans. Our findings provide new morphological support for the hypothesis that Trilobita belongs to the mandibulate stem lineage and bring new insights into the ancestral morphology and evolutionary history of the euarthropod central complex. More broadly, our findings offer direct anatomical evidence of the key sensorimotor system innovations that originated during the Cambrian Information Revolution.

Paraspinal muscles support and move the spine, and reduced size and density of these muscles is associated with increased risk of vertebral fracture and back pain. We test the hypothesis that subsistence-oriented forager-horticulturalists (Bolivian Tsimane) have larger and denser paraspinal muscles, and smaller age-related differences in paraspinal muscle size and density, than more sedentary post-industrialized suburbanites (Framingham Heart Study [FHS]). We measured cross-sectional area (CSA) and density of trunk extensors and psoas major at L3 from computed tomography scans of 411 FHS and 330 Tsimane (50% women, aged 50-90 years). Sex-stratified linear regressions determined associations of muscle CSA and density with age and cohort (Tsimane or FHS [ref.]), with and without an age-cohort interaction term. Covariates were height and weight. After adjusting for covariates, Tsimane women had 11-14% larger and denser extensors, and 17-26% smaller and less-dense psoas major, than FHS women. Tsimane men had 8-21% smaller and less-dense psoas major than FHS men. The age-related difference in extensors density was 41% smaller for Tsimane men than for FHS men. There were no other significant age-related cohort differences in muscle CSA or density. We observed a different balance in the size and density of trunk extensors versus psoas major between Tsimane and FHS. A possible explanation for this observation is cohort differences in the amount and types of physical activity, although future work is required to determine the influence of lifestyle and environmental factors on paraspinal muscle size and density.

Three-dimensional imaging is rapidly increasing the scale of biological morphometric datasets, but establishing automated anatomical correspondence remains a major bottleneck. Conventional point-cloud registration methods transfer landmarks efficiently, yet poor initialization or high morphological variation can produce geometrically plausible but biologically implausible deformations. Here, we developed a statistical-shape-model (SSM)-regularized coherent point drift (CPD) estimator for automated morphometric landmark transfer. It replaces the unconstrained non-rigid CPD update with optimization in SSM coefficient space and regularizes coefficients by inverse population covariance, so deformations along low-variance biological modes are penalized more strongly. We compared pose-matched conventional CPD, SSM-CPD, and SSM-CPD with refinement in leave-one-out atlas experiments on laboratory mouse skulls (N=249), primate molars (N=115), and squamate vertebrae (N=254) under aligned, mildly perturbed, and severely perturbed initialization. We evaluated these formulations in terms of landmark accuracy, failure-tail behavior, runtime, and preservation of morphometric structure. Under severe perturbation, we also evaluated an integrated recovery procedure that optimizes the starting template, recovers global pose, and then applies SSM-CPD with refinement. Severe misalignment produced widespread failures in pose-matched conventional CPD, whereas SSM regularization reduced median error from 13 - 21% to approximately 1 - 3% of the target bounding-box diagonal. This improvement was highly consistent for skulls and molars. In the morphologically complex vertebral dataset, however, low median error concealed a substantial minority of catastrophic failures (28% of vertebrae above a 6% error threshold), which integrated recovery reduced to 4%. Across datasets, automated landmarks preserved the major axes and pairwise relationships of morphometric shape space, and local CPD refinement improved covariance preservation. Embedding learned biological variation directly into registration broadens the range of starting conditions and anatomies for which automated landmark transfer remains reliable. Our results also show that failure tails and downstream morphometric fidelity must be assessed alongside median accuracy when automated landmarking methods are evaluated.

Lineages that have undergone little apparent evolutionary change over long geological intervals are often described as "living fossils", with iconic examples including the coelacanth, tuatara, and Ginkgo. However, studies of morphological and molecular data have produced mixed evidence regarding whether living fossils have evolved more slowly than their relatives. We tested for reduced evolutionary rates using a sister-pairs approach, analysing a comprehensive set of morphological and molecular data for 45 living-fossil lineages and their close relatives across the Tree of Life. Living-fossil lineages showed substantially reduced rates of morphological evolution, whereas evidence for reduced molecular rates was weaker. To test whether this pattern could be explained by species poverty, we also analysed other depauperate lineages (species-poor clades contrasted with species-rich sister clades) and found a weak tendency towards reduced molecular rates but without a corresponding slowdown in morphological evolution. Our study reveals that the living-fossil concept captures a quantifiable pattern of prolonged morphological conservatism, but also suggests that living fossils lack other shared biological or life-history characteristics that might unite them as a group.

Urbanisation imposes novel environmental conditions that can shape phenotypic traits including plumage, yet, its effects on characteristics beyond colouration remain poorly understood. Feather microstructure plays an important role in thermoregulation, protection from solar radiation and waterproofing, and its expression is sensitive to environmental conditions. Here, we quantified feather barbule density, a key component of feather microstructure, in great tits (Parus major) from five urban-rural paired populations across Europe to assess its association with urbanisation, temperature and elevation. We found no significant effects of the urban environment or temperature. However, feather microstructure significantly varied between the studied areas and was positively associated with elevation and ultraviolet B (UV-B) radiation, suggesting a potential adaptive response to variation in solar radiation.

Social relationships are strongly linked to health and biological fitness in humans and other animals. For example, social support may buffer individuals against harsh ecological conditions. Reciprocally, better health may allow for greater investment in social relationships. We examined links between health and social behaviour in an endangered population of northern bottlenose whales (Hyperoodon ampullatus). Using photo-identifications spanning 37 years, we visually classified whether individuals had poor body condition on each day they were observed (N = 3,086 whale-days). We then used generalized linear mixed models to explore how body condition varied with sex, age, and within-individual deviations in group size or other social metrics derived from Bayesian social networks. Approximately 11% of individuals were classified as being in poor condition at any time, though the prevalence of poor condition varied across raters and years. Visually assessed condition was successfully validated with morphometric measures, though there was only weak evidence that poor condition predicted apparent survival. Contrary to expectations, whales in poor condition were found in larger groups and had more social partners than their average across the study period. These findings run counter to the expectation that health and social connectedness are positively correlated in social mammals.

Although Cinnamomum includes globally important spice and medicinal plants, genome-scale knowledge of Cinnamomum species native to Japan remains limited. The Nansei Islands are an approximately 1,200-km-long archipelago located between the Pacific Ocean and the East China Sea and represent a key biogeographic region shaped by deep-sea barriers such as the Tokara and Kerama Gaps. Here, we integrated chloroplast and nuclear genome data generated in this study with public datasets to examine the phylogenetic relationships of four Cinnamomum species native to the Nansei Islands and ten related species. Chloroplast genomes showed no major structural differences among the four island species, whereas nuclear genome data resolved two distinct evolutionary lineages in the Nansei Islands: C. yabunikkei, C. x durifruticeticola, and C. daphnoides formed a distinct Nansei Islands lineage separate from Sri Lankan, Indonesian, and Taiwanese lineages, whereas C. sieboldii was placed in a Taiwanese lineage. Ancestry estimation indicated that C. x durifruticeticola and C. daphnoides share a common ancestral component, while C. sieboldii exhibited a single homogeneous ancestral component, suggesting reduced ancestral diversity after divergence from a Taiwanese lineage. These findings reveal previously unrecognized genomic differentiation among Japanese Cinnamomum species and suggest that deep-sea barriers, dispersal, and long-term isolation shaped their diversification.

Early high disparity, where clades reach maximum morphological diversity early in their history, is often observed in the fossil record. Since there are currently no suitable methods to model these Early High Disparity (EHD) processes accurately, we introduce a new phylogenetic comparative method (PCM) to address this gap by combining two evolutionary modes--Ornstein-Uhlenbeck (OU) and Early Burst (EB)--leading to an initial surge in disparity until an early maximum, driven by the EB component, after which disparity declines due to OU constraints. We tested the EHD model on ichthyosaur skull length data and found it to be the best fit in nearly all cases compared to other modes. The median time of peak disparity (TEHD) is about 231.1 million years ago, during the Triassic, roughly 12.4% into the ichthyosaur lineage. To explore how environmental and geological factors might influence peak disparity, we performed changepoint analyses on sea level, sea surface temperature, rock area, and sedimentary sequence ages. While changepoints in temperature at 237.2 Ma and sequence boundaries at 238 Ma were closest to TEHD, no significant changepoints appeared to match TEHD.

Kashmir voles (Hyperacrius) are poorly studied fossorial arvicolines of the southern Himalayan slopes, with two species recognized. We sequenced Cytb from voles collected in northern Pakistan and combined these data with existing GenBank sequences, yielding 194 tips spanning nearly all arvicoline genera. We analyzed five data sets with maximum likelihood and compared external measurements among species using principal component and linear discriminant analyses. We recovered five well-supported clades of Hyperacrius separated by ~12-14% p-distance, which we treat as distinct species. We recognize clade 1 as H. wynnei, and populations referred to H. w. traubi are genetically indistinguishable from H. w. wynnei. We recognize H. zygomaticus, new rank (clade 2), and tentatively allocate clade 5 to H. brachelix. Clades 3 and 4 likely represent undescribed species. Morphometric analyses of external measurements clearly distinguish H. wynnei but not clades 3 and 4 (overall classification accuracy 70.3%). Phylogenetic analyses including nuclear genes placed Hyperacrius as sister to a clade containing Microtini, Arvicola, Lagurini, Dinaromys, and Ellobiusini, whereas mitochondrial analyses weakly allied it with Ellobiusini. Its traditional placement in Clethrionomyini is thus unsupported, although its exact position remains unresolved. Alticola is nonmonophyletic, with A. macrotis sequences nested within Clethrionomys, likely reflecting mitochondrial capture.

During biogeographical ancestry (BGA) analysis, an unknown individual for whom ancestry informative markers (AIMs) were analyzed is often assigned to one of a few reference populations. For some individuals, such an assignment is not meaningful: the two parents of the individual may come from different populations, or the individual may come from a population not included in the reference database. In this study, we examined classification methods in combination with 104 AIMs from the VISAGE Enhanced Tool to provide one of the following three possible outcomes: (i) in-reference (IR), with an assignment to one of the reference populations; (ii) admixed (AD), with parents from two in-reference populations; and (iii) out-of-reference (OR), where the individual's ancestry is not sufficiently represented by the reference populations. We therefore compared seven classifiers in combination with class-specific acceptance thresholds for the OR decision. The best-performing method according to our selection criterion was a naive Bayes classifier over pairs of populations, which we call the recent naive Bayes (RNB) classifier. Populations held out of the reference panel were correctly rejected in $84.9\%$ of cases. We further tested the method on $800$ individuals from seven additional groups not included in the reference panel. Only $29.4\%$ of individuals from four non-admixed populations were rejected, resembling joint ancestry with reference populations. In addition, we investigated the prediction accuracy in a newly collected dataset of 28 individuals with self-reported ancestry by analyzing AIMs using massively parallel sequencing.

Madagascar's extinct hippopotamuses were a distinctive component of the island's recently lost megafauna, yet their origins, diversification history, and taxonomy have remained unresolved. Here, we recover the first ancient mitochondrial genomes from Malagasy hippopotamuses, obtaining phylogenetic data from 37 subfossil specimens across seven localities and all three named species, and integrate these data with morphology and total-evidence phylogenetic analyses. The 35 endemic specimens form a single clade most closely related to Hippopotamus amphibius, supporting a single colonization of Madagascar. Within this radiation, five deeply divergent and geographically structured mitochondrial lineages diversified across most regions of Madagascar during the Pleistocene. Tip-dating estimates a remarkably young crown age of ~1.41 Ma, following divergence from mainland Hippopotamus at ~3.24 Ma. These lineages do not map cleanly onto traditional taxonomy: material assigned to H. lemerleispans several divergent clades, while H. laloumena and H. madagascariensis are nested within this broader diversity. Morphology supports some lineage assignments but also points to previously unrecognized diversity, while ancestral-state reconstruction supports early body-size reduction followed by a secondary increase in H. laloumena. Malagasy hippopotamuses therefore represent an exceptionally recent radiation of large-bodied island mammals (and the youngest in Madagascar), showing that Madagascar continued to generate endemic mammalian diversity well into the Pleistocene. Their late-Holocene extinction erased the product of more than a million years of rapid insular diversification, the drivers of which remain enigmatic.

Selective attention, the ability to focus on relevant information while ignoring task-irrelevant information, is crucial for learning. Selective attention has largely been studied in motionless individuals. However, body movements are intrinsic to cognition and particularly frequent in children. To test whether movements explain successful attention control throughout childhood, we tracked head rotations while younger children (6-8 years, n = 25), older children (9-10 years, n = 20) and adults (18-30 years; n = 21) identified a spatial target sound under distraction. Distractor reports decreased with age, confirming a developmental increase in attention control. Importantly, a dissociation of two distinct types of movement emerged: Task-dependent head rotation to lateral targets was similar across age groups. Contrary, task-independent movements before and early during a trial decreased with age and predicted attention accuracy. Thus, movements unrelated to the attention goal do not reflect noise but signify attentional control across development.

Life experiences in fathers can influence phenotypes in the offspring. While it is known that acquired paternal information can be transferred to descendants via the germline, how such transfer operates between sperm and embryo remains unclear. We examined the effects of low protein diet (LPD) on males of different genotypes and their hybrid offspring generated by reciprocal outcross and on their genome. We found that LPD induces genotype-specific effects in fathers and different phenotypes in their male and female offspring. Allele-specific analyses revealed that paternal LPD modifies paternal and maternal transcriptomes differently in early two-cell embryos. Changes in the paternal transcriptome overlap with the sperm transcriptome and are associated with sperm chromatin features, suggesting that sperm-inherited transcripts and chromatin states may shape the paternal transcriptome in the embryo. The maternal transcriptome, by contrast, is linked to sperm-derived small RNAs, consistent with a possible role of sperm-derived small RNAs in modulating maternal transcript levels in the embryo. Analyses of undifferentiated spermatogonia populations suggest that LPD-induced sperm states can originate from early stages of germ cells. These findings provide evidence that paternal experiences can influence the early embryo via processes involving each parental allele differently.

Hard cornification is the developmental process that produces the hardened epidermal structures that characterize tetrapod diversity, including scales, feathers, hair, and claws. This process relies on keratin proteins and additional reinforcing proteins that assemble, cross-link, and strengthen tissues to become mechanically resilient. Although the molecular composition of hard cornified tissues is known to vary among tetrapods, the program is least understood in amphibians. Within frogs, hard cornification occurs in the mouthparts of most tadpoles and in the nuptial pads, hindfoot claws, and foot tubercles of some adults. We characterized transcriptomic profiles of Bombina orientalis tadpoles at three developmental stages before, during, and after mouthpart cornification. We identified 25 keratin genes orthologous to amniote hard keratins in the B. bombina genome, 18 of which are upregulated in cornifying larval mouthparts. Furthermore, two of these hard keratin paralogs and an amphibian epidermal differentiation gene exhibit distinct spatial expression patterns within the jaw sheaths and keratodonts evidenced by in situ hybridization, suggesting they have discrete functional roles. The newly discovered diversity of the hard keratin repertoire in Bombina , and the concurrent expression of many of these paralogs in tadpole mouthparts, suggests that anuran hard cornification is much more complex than adult cornified traits thus far characterized in Xenopus frogs. Our results also challenge the notion that the diversification of keratins in Amphibia is primarily associated with the transition to more terrestrial habits, and may instead be shaped, in part, by the novel feeding apparatus of aquatic tadpoles.

Regulatory variation is widely recognized as a driver of evolution within and between species, including in the primate lineage. Evidence for potentially adaptive regulatory evolution in primates comes from both sequence-based tests for selection and cross-species comparisons of gene expression. However, comparative studies of gene expression have tended to focus on baseline cellular states, despite growing evidence that context-dependent gene regulation is important. Here, we investigate the evolution of gene expression in lymphoblastoid cell lines (LCLs) from 6 primate lineages (humans, chimpanzees, bonobos, western gorillas, Sumatran orangutans, and rhesus macaques) using data from both control and 5 perturbed cellular conditions (n=22-25 individuals per condition). Analyzing baseline expression data, we find that genes involved in immunity often exhibit lineage-specific transcriptional shifts, and that these genes occur near regions with lineage-specific sequence acceleration and lineage-specific epigenomic states. However, the subset of genes that respond to in vitro perturbations have evolved at significantly higher rates than environmentally-insensitive genes, a finding we replicate in an independent four-species dataset of primary white blood cells (n=26). Together, our findings support the argument that focusing on baseline states alone provides a biased view of gene regulatory evolution, and suggest that studying gene regulation in dynamic contexts can enrich our understanding of primate evolution.

Compositional syntax - the ability to combine meaningful units into structured sequences whose overall meaning depends on both the meanings of the parts and their combination - is a hallmark of human language. Recent work suggests that chimpanzees may exhibit similar capacities, raising the possibility that compositionality predates the emergence of language in humans. However, a key challenge has been to exclude simpler explanations in which call combinations are processed as independent utterances rather than integrated structures. Here, we address this issue using a source-manipulation paradigm in playback experiments with captive chimpanzees (n = 6) at the Sigean African Reserve, France. We presented sequences combining alarm and recruitment calls either from a single source (one loudspeaker) or from two spatially distinct sources (two loudspeakers). If sequences are interpreted as independent calls, responses should not differ between conditions. Instead, chimpanzees responded significantly faster when sequences originated from a single source. Bayesian survival analyses confirmed this effect. These results support a non-trivial form of compositionality, whereby the meaning of call combinations depends on their integration into a unified sequence produced by a single individual. Our findings provide experimental evidence that key features of linguistic compositionality may have evolutionary roots in the primate lineage.

Skeletal dysplasias can occur in animals, including humans, and are attributed to both environmental and genetic causes. As our understanding of the distinguishing phenotypes of different skeletal dysplasias improves, alongside elucidation of the mechanisms and genetic networks that contribute to skeletal development, identifying and controlling genetic mutations in livestock populations has become increasingly possible, including approaches beyond selective breeding alone. Within Scotland an 'achondroplasia' syndrome that was first described in the South Country Cheviot breed of sheep (OMIA 000004-9940), has continued to occur periodically within the same geographical locations, despite culling of affected lambs over many decades. Affected lambs exhibit rounded heads, short legs, a loss of hooves, abnormal stance and short twisted tails. As the condition has been observed in flocks and as progeny from specific ram and ewe matings, and cases have been co-twinned to normal siblings, this suggests a syndrome caused by a recessive genetic condition rather than environmental factors. Through homozygosity mapping and whole genome sequence analyses, we identified a 10 bp deletion in the open reading frame of the PRICKLE1 gene that encodes prickle planar cell polarity protein 1 and is essential for normal vertebrate embryonic development. We identified additional phenotypes present in affected lambs, including a loss of phalanges, neural tube defects, vertebral and articular chondrocyte defects, hair orientation differences and polycystic kidneys, consistent with a loss of PRICKLE1 function based on comparisons to mouse mutant lines. These findings enabled the development of a simple DNA test for the identified PRICKLE1 mutation that was made commercially available to facilitate the elimination of this genetic disease from sheep flocks.